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Abstract 



By simulating a nonequilibrium coupled map lattice that undergoes an 
Ising-like phase transition, we show that the Lyapunov spectrum and related 
dynamical quantities such as the dimension correlation length are insensi- 
tive to the onset of long-range ferromagnetic order. As a function of lattice 
coupling constant g and for certain lattice maps, the Lyapunov dimension 
density and other dynamical order parameters go through a minimum. The 
occurrence of this minimum as a function of g depends on the number of near- 
est neighbors of a lattice point but not on the lattice symmetry, on the lattice 
dimensionality or on the position of the Ising-like transition. In one-space 
dimension, the spatial correlation length associated with magnitude fluctua- 
tions and the length are approximately equal, with both varying linearly 

with the radius of the lattice coupling. 
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I. INTRODUCTION 



Laboratory experiments |jT| and numerical simulations can now systematically explore 
sustained homogeneous nonequilibrium systems of quite large aspect ratios (T < 1000) 
which possibly approximate a thermodynamic limit of infinite system size. These advances 
raise the theoretical question of identifying order parameters for analyzing and classifying 
spatiotemporal chaotic states so that a quantitative comparison can be made between theory 
and experiment ||^ . The most appropriate order parameter for a given nonequilibrium system 
is presently not well understood although numerous possibilities have been studied. Some 
order parameters, such as spatial correlation lengths obtained from exponentially decaying 
correlation functions, emphasize the average spatial disorder and have been widely used in 
condensed matter physics 0]. Others such as the metric entropy and the Lyapunov fractal 
dimension P] are familiar from nonlinear dynamics and emphasize the average temporal 
disorder or dynamical complexity arising from the geometric structure of an attractor in 
phase space. 

We would then like to know whether these different kinds of order parameters are related 
and whether there is a need for new order parameters. As an example, does knowledge of an 
easily measured correlation length give information about the fractal dimension, which is dif- 
ficult to estimate from experimental time series [^? That a relation between spatial disorder 
and dynamical complexity may exist is suggested by the prominence in many nonequilibrium 
states of defects whose dynamics often determine the average spatial disorder . Exam- 
ples of defects are amplitude holes in one-dimensional complex Ginzburg-Landau equation 
(abbreviated below as CGLE) |9Hr^, domain walls and droplets of opposite spin in coupled 
map lattices (CMLs) with an Ising-like transition |jl3|,|l4|, vortices in the two-dimensional 
CGLE IQ], and spirals and centers in the recently discovered spiral-defect chaos state in 
thermal convection [0. The nucleation, motion, and annihilation of defects are important 
features of the chaotic dynamics and so dynamical quantities such as the fractal dimension 
may be related to their spatial statistics |]I6|-[T8|]. Complicating this simple picture is the 
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fact that not all fluctuations are associated with defect motion, e.g., phase fluctuations in 



the CGLE |TT| , p!9| . The fractal dimension may then be larger than that suggested by defect 
statistics. 



In recent work p7|JT^JT9[] , Egolf and Greenside have explored the relation between tem- 



poral and spatial disorder for the CGLE on a large periodic interval W^. Observing that 
sufficiently large chaotic systems become extensive so that a fractal dimension D grows 
linearly with volume size V oc L'^ (where L is the system size and d is the spatial dimen- 
sionality) [l7|,p0 H2^JT7|j23|j2^ , they calculated the dimension correlation length which is 



defined [0] in terms of the intensive dimension density 

6 = lim D/V, (1) 

L— >oo 

by the equation 

= r^/^. (2) 

The length C,s can be interpreted crudely as a characteristic size of dynamically indepen- 
dent subsystems or a characteristic range of chaotic fiuctuations. Near a transition from 
phase- to defect-turbulent states PT| , P5| , Egolf and Greenside found that the length was 
approximately equal to, and had the same parametric dependence as, the spatial correlation 
length given by the magnitude fiuctuations of the Ginzburg-Landau field ||12[. Over the 
same parameter range, the correlation length arising from phase fiuctuations (and also 
of the Ginzburg-Landau field itself) was found to increase to quite large values, suggesting 
that the chaos fiuctuations measured by were short-ranged and decoupled from the phase. 

In this paper, the one- dimensional investigations of Egolf and Greenside ||r?| , p!2| are ex- 
tended by examining similar questions of how spatial disorder and dynamical complexity 
are related for spatial dimensionality d = 2 and d = 3. We study a class of dissipative 
coupled map lattices that undergoes an Ising-like ferromagnetic-ordering transition as 
a lattice coupling constant g (defined below in Eq. (^) is increased through a range of finite 
positive values, corresponding to a transition from a high-temperature paramagnetic phase 



to a low-temperature ferromagnetic phase. Instead of the point-hke space-time defects of 



the one-dimensional CGLE |[T^, long-lived defects occur in the form of domain walls and 
droplets involving regions of opposite sign. At a critical transition point g = Qc, the usual 
two-point correlation length ^2 diverges to infinity while the magnetization (the average of 
the signs of all lattice variables) bifurcates from a zero to nonzero value. The main attraction 
of the Miller-Huse model fl^ is that a transition with a diverging correlation length occurs 



from one chaotic state to another. This permits a careful comparison of different length 
scales near the transition point. 

By calculating the Lyapunov spectrum and related dynamical quantities such as the 
dimension correlation length ^5, we show that correlations in chaotic fluctuations near the 
ordering transition have a short range and are decoupled from the diverging long-range order 



measured by the length ^2- As was the case for the Id periodic CGLE [0, the length 
and the correlation length arising from magnitude fluctuations of the CML variables 
both turn out to be short — here about one lattice spacing — but a quantitative relation can 
not be determined since these lengths do not vary substantially with parameters. In one- 
space dimension, a more substantial variation in these quantities is obtained by increasing 
the radius r of the coupling from nearest to rth-nearest lattice neighbors. With increasing 
radius r, the lengths C,s and then increase approximately linearly and with slopes related 
by a factor of order one. This suggests that the two lengths may be related and that r is 
important in determining the length scale of dynamical fluctuations. 

Several dynamical quantities such as the largest Lyapunov exponent Ai, the metric en- 
tropy density h = limj;^^oo H/V, and the Lyapunov fractal dimension density 6 attain mini- 
mum values near — but distinctly not at — the critical transition point g = Qc- At first glance, 
minima in quantities such as h or 6 seem counterintuitive since the onset of ferromagnetic 
ordering should correspond to increased correlations, i.e., decreased dynamical complexity 
and a decreased dimension density or entropy density. Such a monotonic decrease in the 
metric entropy is, in fact, observed for a non-dissipative equilibrium CML that undergoes an 
Ising-like transition . But because the dimension correlation length ^,5 is quite short in the 



Miller-Huse model, we argue below that dynamical quantities are only sensitive to nearest 
neighbor dynamics. As the coupling constant g increases, the discrete Laplacian eventually 
becomes antidiffusive, magnifying rather than reducing short-wavelength structure, and the 
dimension density and entropy density start to increase. 

The role of coupling lattice neighbors is demonstrated with calculations on periodic 
2d hexagonal and 3d cubic lattices. For these cases, is again about one lattice spacing in 
size and the extrema in dynamical quantities occur at a value g \/{n + l) determined by 
the number n of nearest neighbors (n = 4 for the square lattice, n = 6 for the hexagonal 
and cubic lattices). The positions of the extrema do not coincide with the bifurcation of the 
magnetization and do not seem to depend on the lattice symmetry and dimensionality. The 
origin of these minima remains to be explained. 

These results, together with previous work on the Id CGLE and with some unpublished 



work on CMLs with algebraic decay of spatial correlations p4[, suggest that the dimension 
correlation length will typically be short so that chaotic fluctuations are decoupled from 
long-range spatial order as measured by correlation functions. This leads to the qualitative 
conclusion that defects, although a striking visual feature of spatiotemporal chaos, are not 
the source of complexity that leads to large fractal dimensions and to small dimension corre- 
lation lengths. The physical meaning and utility of the length remains to be understood 
and further studies on different kinds of systems will be useful. 

The rest of this paper is organized as follows. In Section we define the coupled map 
lattice and discuss some details of how its Lyapunov spectrum is calculated using a CM-5 



parallel computer p6| . In Section |T|, we discuss various results of our simulations, especially 
the dependence of order parameters on the lattice coupling constant and on the symmetry 
and dimensionality of the lattice. Finally, in Section |V|, we summarize our results and relate 
them to other recent research. 
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II. METHODS 



In this section, we define the mathematical models used in our simulations and discuss 
some details about how the Lyapunov spectrum and spatial correlation lengths were calcu- 
lated numerically. Since a CM-5 parallel computer played an important role in our being 
able to explore large space-time regions for many parameter values, we also discuss some 
details of how the algorithms were adapted for parallel computation. 



A. Models 

As easily simulated models of spatiotemporal chaos, we consider homogeneous coupled 
map lattices (CMLs) in which the same chaotic map (p{y) is associated with each point of a 
finite periodic lattice and for which nearest neighbor maps are coupled linearly by diffusion. 
CMLs have a significant advantage over partial differential equations of being analytically 



amenable [27| and easier to simulate on a computer. CMLs have the drawback that their 
solutions can not generally be related quantitatively to experiment and they may not have 
universal critical properties near transitions p8[. 



We study CMLs suggested by recent work of Miller and Huse who analyzed the long- 



wavelength properties of a two-dimensional CML that orders ferromagnetically, in analogy 
to the equilibrium Ising model Following these authors, we use a lattice map (j){y) = 



—(f){—y) with odd symmetry so that domains of opposite "spin" or sign arise. The odd 
symmetry is a necessary but not sufficient condition [^] for obtaining an Ising-like transition 
in which the magnetization (defined below in Eq. (H)) bifurcates to a nonzero value as a 
coupling constant g is varied. Following Miller and Huse, we choose (j){y) to be a piecewise- 
linear map with slope of constant magnitude greater than one, for which only chaotic states 
of constant measure exist in the absence of lattice coupling. 

If yj denotes the variable at spatial site i at integer time t (with t = 0, 1, ■ ■ ■), then the 
rule for updating each lattice variable to time t + 1 is given by fl^ : 
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yl^' = m+9T.{^iy',)-<f^iyl))^ (3) 

Hi) 

where the parameter g is the spatial couphng constant. The sum goes over indices 
that denote nearest neighbors sites of site i, e.g., the four nearest neighbors on a square 
two-dimensional lattice or the six nearest neighbors in a 2d hexagonal lattice or 3d cubic 



lattice. For most of our calculations, we used the same local map as in Ref. [|14 



-2 - 3y for -1 < y < -1/3, 
(p{y) = { 3y for -1/3 < y < 1/3, (4) 
2-3?/ for 1/3 < y < 1, 

with a slope of constant magnitude equal to 3. To understand the competition between local 
chaos and diffusion (which decrease and increase spatial correlations respectively), we also 
used a more weakly chaotic map with a slope of smaller constant magnitude equal to 1.1: 

-l.lx - 1.0 for X < 0, 

<P{y) = { ~ (5) 

I -l.lx + 1.0 for x>0. 

Once a local map (j){y) has been chosen, a numerical simulation is specified by the di- 
mensionality of the lattice d, the symmetry of the lattice (e.g., square, hexagonal, or cubic), 
the size of the lattice L (number of sites along an edge), the coupling constant g, the initial 
condition y^, and the total integration time T. In this paper, we used an integer lattice in Id, 
square and hexagonal lattices in 2d and a cubic lattice in 3d. Initial conditions consisted 
of assigning a random uniformly-distributed number in the interval [—0.1,0.1] to each site; 
results were not dependent on the choice of initial conditions provided the integration time 
was sufficiently long. Typical integration times were T = 50, 000 iterations for calculations 
of Lyapunov exponents and T = 150, 000 iterations for calculations of correlation functions. 
We made a few runs with longer integration times of T = 60, 000 and T = 500, 000 to 
check the convergence of the Lyapunov exponents and the correlation functions, respec- 
tively. Different lattice sizes were used depending on which order parameters were being 
studied, typically L < 32 for calculating Lyapunov exponents (which are quite costly to 



compute) and L < 1024 for estimating correlation lengths. Many runs were repeated using 
several lattice sizes to verify the absence of significant finite-size effects. 

For calculations of statistical averages such as the magnetization and correlation func- 
tions, the statistics could often be improved substantially by averaging the results of an 
ensemble of N runs each of duration T, with each run differing only in the choice of initial 
conditions. Calculations indicate that this ensemble average is ergodically equivalent to a 
single integration of duration NT |ll9| , |3l| . For most of the results reported below, we used an 



ensemble average over = 64 runs which could be executed simultaneously and in parallel 
on the vector units of a 16-node partition of a CM-5 computer. 

B. Lyapunov Exponents and Associated Dynamical Quantities 

Some dynamical order parameters can be constructed by combining in different ways 
the Lyapunov exponents Xi associated with a given attractor 0]. For the CML given by 
Eq. with a total of N = L'^ lattice sites, there are N real-valued Lyapunov exponents Aj 
(labeled in decreasing order Ai > A2 > • ■ ■ > Ajv) that characterize the long-time average- 
rate-of-separation of nearby orbits in phase space. From the Aj, we can calculate a Lyapunov 
dimension D given by the Kaplan- Yorke formula [^: 

D = K + -^j2^i^ (6) 

I'^A'+ll i=l 

and calculate an entropy defined by the sum of the positive exponents 0: 

if = E A. (7) 

Ai>0 

The number K in Eq. (^) is the largest integer such that the sum X^^i K of the first K 
exponents is nonnegative; this sum is positive for i^' = 1 if an orbit is chaotic (Ai > 0) and 
is negative for K = N if the dynamics is dissipative and so the sum typically crosses zero 
at an intermediate index K D for a chaotic dissipative system. 

The exponents Aj were calculated numerically by a now-standard numerical method P2 



in which K < N linearizations of the equations of motion, Eq. (|^), are evolved in time. This 



allows one to follow K Lyapunov vectors in a tangent space from which local stretching 
information and the Lyapunov exponents can be extracted. Together with a particular 
nonlinear orbit defined by the equations of motion, a total of + 1 CMLs was evolved to 
calculate K Lyapunov exponents. 

Repeated orthonormalizations of Lyapunov vectors at time intervals T„ are needed to 
prevent floating-point overflow from the exponentially growing values and to prevent inac- 
curacies arising from the loss of linear independence as they fold up along the direction of 
the fasting growing exponent |^2|. For the maps Eqs. @) and (1), we found empirically that 



values Tn ^ 30 gave reasonable results for all lattices studied with the largest value of T„ 
depending on the parameters g and L. Smaller renormalization times did not change the 
values of the Lyapunov spectrum (although the code was more costly to run) while larger 
values led to serious errors due to linear dependence of the Lyapunov vectors. 

The orthonormalizations of tangent vectors consumed most of the computing time on a 
Thinking Machines CM-5 parallel computer. The orthonormalizations require substantial 
communication between different processors since each processor evolves independently only 
a few of the K linearized equations in its own local memory; this communication decreases 
the efficiency of the code. 

Although the communication inherent in the orthonormalization procedure could not be 
avoided, in all other portions of the code communication between nodes was eliminated by 
iterating redundantly an identical copy of the nonlinear CML Eq. (0) with identical initial 
conditions on each processor. In this way, information about the underlying orbit (needed 
when iterating the linearized CMLs) did not have to be communicated from one processor 
to all others at each time step. 

By monitoring the Lyapunov exponents and the Lyapunov dimension as a function of 
time t, we found empirically that an integration time T ^ 50, 000 time steps gave an 
acceptable relative accuracy of better than one percent for calculating the dimension D and 
entropy H for all lattices studied (L < 32). Fig. [I| shows how the dimension D converges over 
time for lattice size L = 16 and for g = 0.202. The dimension curve is noisy with fluctuations 



that diminish slowly over time (the Xi, not shown, have substantially noiser time series). 
Goldhirsch et al have argued that the amplitude of the exponent fluctuations should 
decay as 1/T where T is the total integration time and so one could flt and extrapolate to 
get a better estimate |Q . Extrapolation was not needed in plots like Fig. |1| which already 
give an adequate relative accuracy of better than one percent. 

By repeating plots such as Fig. || for different system sizes L with all other parameters 
held flxed, we found extensive scaling of the dimension D with the volume of the system N = 

for a wide range of parameter values g, an example of which is given in Fig. |^. Fig. |^(a) 
shows that the dimension D increases linearly and extensively with N beyond a system size 
of about L = 9. The Lyapunov dimension density 6 could then be obtained from the slope of 
a least-squares fltted line in the extensive region. The intercept of a least-squares-fltted line 
through the four right-most points was 0.007 which is quite small (and is also approximately 
zero for extensively chaotic solutions of the Id CGLE [1^). There is then the possibility 



that a single dimension calculation for a sufficiently large system may suffice to estimate 
its dimension density. By comparing the dimension per volume D/L"^ with the dimension 
density 6, Fig. ^(b) shows how the extensive regime is approached rapidly and achieved for 
fairly small system sizes L > 9. 

We finish this subsection with two comments about the meaning of the Lyapunov di- 
mension D and about why we chose to calculate D from Lyapunov vectors rather than from 
time series measurements. The reader should recall that there is an infinity of fractal dimen- 
sions Dg (often called the Renyi dimensions) associated with a strange attractor with the 
parameter q varying over the real numbers 0. Since each dimension Dg will be extensive for 
a homogeneous extensively chaotic system, the particular values only refiect the system size 
and are not interesting themselves. Instead, one should define a continuum of corresponding 
intensive dimension densities 6g = lim^^^oo Dg/V to provide a partial characterization of 
such systems. 

For large fractal dimensions (D ^ 5), present computers and algorithms only allow the 
calculation of the Lyapunov fractal dimension Eq. (P) and its corresponding density. The 
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Lyapunov dimension is conjectured to be the same as the Renyi dimension with q = 1 (i.e., 
the information dimension Di 0) and it is not known to what extent the corresponding 
density 6i characterizes the unknown function of densities 6q, e.g., whether it is close to 
the mean value of the 6q. There are one-dimensional maps for which the ratio of Di to D2 
(the two most commonly calculated fractal dimensions) can be arbitrarily large and 
so the variation of the function Sq around its mean value can be large. A perhaps even 
more important question is whether the different dimension densities 5q each have a similar 
dependence on model parameters, e.g., all increasing or decreasing together. Until this issue 
is resolved, the dimension density 61 and the corresponding length = ^^"^ need to be 
interpreted with caution. 

The dynamical quantities D and H are calculated in terms of the Lyapunov spectrum Aj, 
and not in terms of time series yj at a given lattice site i, because of the impractical com- 



putational demands of time series algorithms [36|. While the computational complexity of 



the method based on the Lyapunov spectrum scales algebraically with system volume V 
or dimension D [^, it is well known that the computational complexity of classical time 
series algorithms such as that proposed by Procaccia and Grassberger grows exponentially 
with D, imposing severe restrictions on the largest dimension that can be estimated from 
experimental data. Although it remains controversial what is the practical upper bound for 
clean time series of less than a million points (there are claims from 6 to 20 @]), existing 
time series methods can not treat extensively chaotic systems whose fractal dimensions may 
be in the hundreds (see Fig. 0). 

C. Magnetization and Correlation Lengths 

In addition to the dynamical quantities described in the previous section, we quantify 
the evolution of the CMLs by a "magnetization" M and by length scales measured from two- 
point and mutual information correlation functions. We describe these briefly to indicate 
the method and errors involved. 
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Following Miller and Huse [0, the average magnetization M of the CML is defined by 



a space-time and ensemble average of the signs ±1 of the lattice values 

M = (sign(7/*)) = ^ sign(|/*(p)), (8) 

where the index p labels a particular CML running on processor p. An average over 64 inde- 
pendent CMLs running on a 16-node partition of a CM-5 was typically used. For the local 
maps Eq. @) and Eq. (|^), M undergoes a pitchfork bifurcation from a zero to finite value 
as the lattice coupling g is increased from small values. The critical value Qc at which M 
bifurcates to a nonzero value is unchanged if the values of the site variables y\ are used 
instead of their signs in Eq. (|^). The bifurcation of M to a nonzero value defines the onset 
of ferromagnetic order at g = Qc [|14| as illustrated in Fig. |^. 

To characterize the average spatial disorder, we examined two of many possible defini- 
tions of spatial correlation lengths, one from an exponentially decaying two-point correlation 



function, another from an exponentially decaying mutual information function |3^. (Some 
other correlation lengths are discussed on pages 945-947 of Ref. [0.) The two-point correla- 
tion function was defined in the usual way: 



C,i\X,-X,\) = {{yl-{y)){yl,-{y))), (9) 

where Xj denotes the position of lattice point i and where the brackets (■) denote the 
averaging process of Eq. (1). Given the periodicity of the lattice and the availability of 
efficient parallel Fast Fourier Transforms (FFTs) on the CM-5, we calculated Eq. via 
the Wiener-Khintchin theorem , first obtaining the time-averaged magnitude squared of 



the Fourier coefficients, from which Eq. @ was obtained by an inverse FFT. In most cases, 
there was a substantial region of exponential decay from which the two-point correlation 
length ^2 was obtained by a least-squares fit of the form aexp(— a;/^2); a representative plot 
is given in Fig. ^. The correlation functions and corresponding values of .^2 do not change if 
Eq. is defined in terms of the sign of the variables, sign{yj). 

The two-point correlation functions decay more rapidly as shown in Fig. ^ if the signs of 
the field values yj in Eq. (||) are replaced by their magnitudes \yj\. The correlation length ^^'^^ 

12 



obtained from the initial rapid exponential decay is approximately one lattice spacing and 
changes little when the coupling constant g is varied over a large range, including near the 
bifurcation point g = gc- For dimensionality = 1, we show below that this short length 
scale ^^'^^ is related to the dimension correlation length and that both vary linearly with 
the radius r of neighboring lattice sites that are coupled together spatially. 

As a possible alternative for characterizing the spatial disorder of a nonlinear system, 
we also calculated a correlation length C,i based on the exponential decay of the mutual 
information function /(AX) |^8| of the variables yj. Fig. ^ shows the exponential decay of a 
mutual information function /(AX) for two-dimensional square lattice of size L = 256 and 
for the parameter g = 0.204. Again if the magnitudes of the field variables are used when 
calculating /(AX), the exponential decay is much more rapid, with ^^^^ being approximately 
one lattice spacing for a wide range in g. 

Although there is not yet a compelling theoretical reason to prefer ^/ over other corre- 
lation lengths such as ^2 105 former has the distinction of depending nonlinearly on the 
dynamics and so may depend on details that are missed by ^2- For this reason, an increasing 
number of scientists have reported correlation lengths in terms of ^/ pn| , p5| , ^ . As shown 
in Fig. 1^ for the two-dimensional Miller- Huse CML on a square lattice, the length scales ^2 
and ^/ are linearly related over a substantial dynamical range near the ferromagnetic tran- 
sition. At least for the present models, these lengths are equivalent measures of spatial 
disorder and we report values only for ^2 below. 



III. RESULTS AND DISCUSSION 

In this section, we discuss our calculations of the Lyapunov spectrum and of correlation 
lengths. Our goal is to explore how spatial disorder (as characterized by the two-point 
correlation length or by the mutual information correlation length) is related to dynamical 
complexity (as measured by the intensive dimension density and by the dimension correlation 
length Eq. (^) and to investigate how these order parameters vary near the nonequilibrium 
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transition point g = Qc a.t which the magnetization bifurcates to nonzero values (Fig. |^). 

Results for the 2d square lattice are given first, followed by results for lattices with 
different symmetries and dimensionalities. We do not address issues related to critical ex- 



ponents of these different models which have been discussed by Miller and Huse [|T^] and 
more recently by Marcq and Chate PS|. Related interesting results on similar CMLs have 
also recently been reported by Boldrighini et al . 



A. Results for Two-Dimensional Square Lattices 

For the two-dimensional periodic square lattice with map Eq. Fig. ^ shows that there 
is a bifurcation at Qc — 0.205. This bifurcation corresponds to the onset of long-range order 
of the lattice variables as demonstrated by the divergence of the two-point correlation 
length ^2 as g approaches Qc (Fig. §(a)). Over this same parameter range, the dimension 
correlation length C,s varies smoothly (Fig. §|(b)), deviating by less than four percent from a 
value of one lattice spacing and attaining a maximum value close to where the correlation 
length diverges. The Lyapunov spectrum of exponents also varies smoothly from one side of 
the transition to the other as shown in Fig. ^. We conclude that chaotic fluctuations have 
a short range, are decoupled from the onset of long-range order measured by and that 
the spectrum of exponents is at most weakly dependent on the onset of long-range spatial 
order. 

To understand further how various dynamical quantities change near the transition point, 
we have plotted in Fig. |10| the variation of the Lyapunov fractal dimension density 6, of the 
metric entropy density h, and of the largest Lyapunov exponent Ai across the ferromagnetic 
transition for a lattice of size L = 16, which is already extensively chaotic according to Fig. |^. 
As was the case for the length in Fig. ^(b), these quantities change by only a small amount 
through the transition (at most by 20%) and all go through a minimum close to, but distinct 
from, the ferromagnetic transition at Qc ~ 0.205. This result was surprising to us since one 
consequence of coupling neighboring maps more strongly (increasing the parameter g) would 
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intuitively be to increase correlations between their dynamics, which should decrease both 5 
and h. Fig. p!0|(a) indicates that roughly one quarter of the maximum number of degrees 
of freedom disappear when the lattice attains its minimum dimension density of 5 ~ 0.746. 
(An upper bound of 5 = 1 is set by the integer lattice spacing.) 

It is not clear why the fractal dimension density and other dynamical quantities have 
extrema near g = 0.20. For an equilibrium non-dissipative CML of Ising dynamics on a two- 
dimensional square lattice, Sakaguchi did not find a local minimum in the entropy H 
but instead found a monotonic decrease consistent with the analytical solution of the spin- 
1/2 Ising model on a square lattice ||29|. One explanation for the extrema may be that the 
dissipative linear coupling in Eq. (|^) becomes antidiffusive for g > 1/5, enhancing rather 
than damping short- wavelength fluctuations and so decorrelating nearby lattice variables. 

The issue is somewhat more subtle than this because the existence of the minimum 
depends also on details of the local map (j){y) in Eq. (^. For the less chaotic lattice map 
Eq. d^) with slope of constant magnitude 1.1, an Ising-like transition still occurs as shown 
by the bifurcation of the magnetization near g^ ~ 0.168 in Fig. |TT]. Fig. ^ now shows 
that the Lyapunov dimension density S and entropy density h decrease monotonically as the 
parameter g is increased, with the largest exponent Ai remaining constant. 

B. Results for Other Lattices 

Figures ^ i(a), and |T0| suggest that the extrema of dynamical quantities may be related 
to the ferromagnetic transition. On the other hand, the short dimension correlation length 
in Fig. Kb) contradicts this by implying that chaotic fluctuations occur over a length scale 
that is short compared to the ferromagnetic ordering. To understand this further, we have 
explored CMLs of different symmetry and dimensionality and found that the near-proximity 
of the extrema with the transition is a coincidence for the two-dimensional lattice with square 
symmetry. More generally, the positions of extrema seem to be determined simply by the 
number of nearest neighbors n, and not by the symmetry or dimensionality of the CML or 
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by the position of the magnetization bifurcation point. 

Fig. ^ summarizes calculations for a two-dimensional periodic hexagonal lattice by plot- 
ting the dependence of magnetization M and of dimension correlation length S,5 on the cou- 
pling constant g. The magnetization bifurcates to a nonzero value at 0.120 which is 
a smaller value than that for the square lattice since the larger number of nearest neigh- 
bors (six versus four) increases the effective strength of the diffusive coupling. The relative 
difference between the transition aX g = gc and the positions of the extrema in and in 
related dynamical quantities is substantially larger than was the case for the 2d square lat- 
tice. For the hexagonal lattice, extrema in quantities like the length ^5 occur at a value close 
to g = 1/ {n + 1) 0.143 where n = 6 is the number of nearest neighbors. 



A similar result is found for the same CML on a 3d cubic lattice, as shown in Fig. [14. 
The transition at g ~ 0.11 occurs at a value close to but smaller than the value on the 
hexagonal lattice. Extrema in the dynamical quantities like again occur at a value close 
to l/(n + 1) with n = Q. 

That the positions of the extrema of dynamical quantities is dependent primarily on the 
number of nearest neighbors is a consequence of the nearest-neighbor diffusive coupling in 
Eq. and of the fact that, for the local map Eq. (^), the chaos is sufficiently strong to make 
the dimension correlation length ^5 quite small, about one lattice length. That the positions 
of the extrema seem to be given quantitatively by the specific formula = l/(ra + 1) is more 
delicate to understand but may be related to the fact that the discrete Laplacian operator 
changes from diffusive to antidiffusive behavior at this value. The value g = l/(n + 1) is 
the value for which the weight of each of the neighbors is equal to the weight of the central 
lattice site to be updated. 

For all CMLs that we studied, the dimension correlation length C,s was about one lattice 
spacing in size and this was also the "radius" of the diffusive coupling in Eq. (^. This 
suggests that the length may be determined by the spatial extent of the diffusive coupling, 
becoming larger as more sites are coupled to a given site. This conjecture was tested in one- 
space dimension by coupling together, with equal weight g, all lattice variables within a 
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radius r of a given site i: 



i+r 



vt"' = Hvl) + 9 E {<Piy',) - <Piyl)) ■ (10) 

j=i-r 

For this one- dimensional periodic CML with the lattice map Eq. (^), the dynamics varies in 
a complicated way with radius r. For most initial conditions, chaos was found for smaller 
radii (r < 6). For larger radii r > 6, the transients lasted much longer and the asymptotic 
dynamics was periodic. As an example, for r = 10 the dimension as a function of time 
initially reached a value D = 50 even after 5000 transient iterations were skipped; however, 
the dimension then decreased steadily to zero over the next 30,000 iterations. We believe 
that this asymptotic periodic behavior is a finite-size effect. For a sufficiently large system 
size L, with the crossover length increasing with the radius r, the asymptotic state should 
be chaotic. 

For dimensionality d = 1, the dimension correlation length varies more strongly with 
increasing radius r than with coupling constant g, which allows several different length scales 
to be compared. Fig. shows that the two-point correlation length ^^'^^ obtained using 



the magnitudes \yj\ of the field values has approximately the same linear dependence on the 
coupling radius r as the dimension correlation length C,5- In addition, these two length scales 
are the same order of magnitude. The two-point correlation length ^2 obtained using the 
actual field values is larger than C,s and does not have the same simple linear dependence on 
r. The situation is then similar to results found for spatiotemporal chaotic solutions of the 



Id periodic CGLE |12| in that the spatial correlation length of fiuctuations in the magnitude 



of a field provides a way of estimating the length ^,5- 



IV. CONCLUSIONS 



In this paper, we have extended recent calculations |0,[T2|JH] concerning the relation 
between spatial disorder and dynamical complexity of a sustained homogeneous nonequi- 
librium system from dimensionality d = 1 to dimensionalities d = 2 and d = 3. This 
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was accomplished by choosing a coupled map lattice, Eq. (Q), that underwent an Ising-like 
transition with diverging two-point correlation length as a parameter g was varied |T^. By 
comparing various length scales such as the two-point correlation length ,^2, the dimension 
correlation length ^5, and the two-point correlation length of magnitude fluctuations 
near the transition point, we were able to show that the lengths ^5 and ^^^^ were short, of 
order one lattice spacing, even as the length ^2 diverged to infinity. In agreement with pre- 



vious work [|T^], the chaotic fluctuations are decoupled from the average long-range spatial 
order. The correlation length of magnitude fluctuations seems to provide an effective way 
to estimate the size of the length ^5. 

Our calculations of the Lyapunov spectrum and related quantities such as the Lyapunov 
dimension density 5, metric entropy density /i, and the largest Lyapunov exponent Ai show 
that the onset of long-range spatial order (diverging ^2) does not affect dynamical invariants, 
which vary smoothly and weakly through the transition point g = gc- Thus the average spa- 
tial disorder (measured by ^2) does not determine dynamical complexity (measured by C,s)- 
Rather surprisingly, the intensive densities S and h go through a minimum near the tran- 
sition point so that the onset of long-range order does not correspond to a decrease in 
complexity. By examining CMLs of different symmetry and of different dimensionality d, 
we showed that the positions g of the extrema were determined by the number n of neigh- 
bors nearest to a given lattice site (with g ^ l/{n + 1)) but not by the symmetry or by d. 
This result can be understood as a consequence of the extremely short dimension correla- 
tion length ^s, about one lattice size, so that lattice variables are independent except when 
they are nearest neighbors. We believe that the minima in 6 and in h occur approximately 
when the discrete Laplacian in Eq. (^ becomes antidiffusive with increasing parameter g. 
Short- wavelength fluctuations are then magnified instead of damped, decreasing correlations 
between neighboring sites. 

Some of our results concerning extrema in dynamical quantities have been independently 



obtained by Boldrighini et al PD| although these authors worked with extensive, rather than 



intensive, quantities and they did not determine whether their calculations corresponded to 



extensively chaotic regimes. Boldrighini et al investigated CMLs of the form Eq. for 
dimensionahties d = 1 and d = 2 but with some new local maps (piy). Besides also finding 
extrema in the Lyapunov fractal dimension, Boldrighini et al showed that the odd symmetry 
of the map (j){y) was not a sufficient condition for the magnetization to bifurcate to a nonzero 
value. Using a strongly chaotic local map with slope of constant magnitude equal to 5, they 
also showed that the magnetization M did not bifurcate to a nonzero value if the local map 



were made sufficiently chaotic compared to the ordering caused by diffusion. In Section |IIIB 
we used the less-chaotic local map Eq. @ with slope of constant magnitude 1.1 to show 
that the dimension density 6 can decrease monotonically without a minimum even when the 
magnetization M bifurcates to a nonzero value. The dependence of these minima on details 
of the local map is not yet understood and should be pursued with further studies. 

The small values of C,s in the present CMLs, in the Id CGLE [Q, and in a nonequilibrium 



CML with algebraic decay of spatial correlations ||2^ have several interesting implications. 
One is that many previous laboratory experiments |1[] and numerical simulations 0] concern- 
ing spatiotemporal chaos are likely already extensive so that it is meaningful to talk about 
dimension and entropy densities (see also the earlier paper by Bohr A second impli- 

cation is that the dimension correlation length C,5 may not be a useful order parameter for 
future studies of spatiotemporal chaos since it depends only weakly on parameters. A third 
implication is that the short value of C,5 suggests the nonexistence of macroscopic chaotic 



states for dynamics with local interactions, a point already made by several researchers 
Finally, we speculate that C,s is the length scale below which one can replace chaotic fluc- 
tuations by a white noise source when trying to develop a hydrodynamic (long wavelength) 
description of spatiotemporal chaos (see the discussion on pages 953-954 in Ref. 0). 

The short values of raise the question of what determines this length scale. Our 
calculations on the Id CML Eq. (|T0|) with a variable radius of coupling r suggest that 
the length is determined partly by the length characterizing magnitude fluctuations 
although the reason for and the generality of this correspondence is not understood ||T^ . The 
length ^5 is also related to the radius r over which nearby lattice sites are coupled together 
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(Fig. |T5]). Further calculations with different kinds local maps and of diffusive operators and 
for different values of r should provide further insight. 

It is appropriate to finish with a discussion about the relevance of these results for 
laboratory experiments. As discussed at the end of Section |11B| , it does not seem possible 
in the near future to calculate the Lyapunov spectrum, the fractal dimension, or the fractal 
dimension density of a high- dimensional extensively chaotic experimental system for which 
only time series measurements are available ^j. Our success in calculating these quantities 
was a result of having explicit knowledge of the dynamical equations which could then be 
integrated numerically on a powerful parallel computer using algorithms whose complexity 
only grew algebraically with the dimension D ||3^. For many laboratory experiments, a 
quantitative mathematical description is either lacking (e.g., for chemical reactions) or, if 
known, is too difficult to work with numerically (e.g., the five three-dimensional Boussinesq 
equations describing buoyancy- induced convection in a large-aspect-ratio container |^). 

Our calculations in Section ^ suggest that one possible way to estimate the dimen- 
sion correlation length C,s may be to calculate the correlation length of some function of 
the physical fields, e.g., the field magnitude. Another possibility will be to discover and 
to validate algorithms that can calculate the intensive dimension density Eq. (|I|) directly 



from time-series measurements that are localized in space in lieu of calculating a large 
extensive fractal dimension D and then dividing by the extensive system volume V. Several 
steps have been taken in this direction but a theoretical foundation has not yet been 
established nor have the numerical algorithms been adequately tested. 
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FIGURES 

FIG. 1. Lyapunov fractal dimension D, Eq. (|6|), versus iteration number t for a periodic 
2d CML on a square lattice with local map Eq. with coupling constant g = 0.202 and with 
lattice size L = 16. Fluctuations in D damp out slowly with increasing time, giving a relative 
accuracy of about 0.1% in the dimension (here estimated to be D = 191.2). 

FIG. 2. Lyapunov fractal dimension D versus the number of lattice points N = L'^ (which 
is also the system volume) for a 2d periodic CML on a square lattice with local map Eq. (g) and 
coupling constant g = 0.199. Each dimension value was obtained from a plot similar to Fig. |l| 
over a time scale of T = 50, 000 iterations, (a) The Lyapunov dimension increases linearly and so 
extensively with for > 100. A least-squares fitted line through the four right-most points gives 
a Lyapunov dimension density (slope) of (5 = 0.746 and an intercept of 0.0077. (b) The normalized 
deviation {D / L? — 5)/ 5 of the dimension density 5 from the dimension density predicted from 
the dimension D per volume Li^ illustrates the rapid onset of extensive chaos for small system 
sizes L > 9. 

FIG. 3. Average magnetization M (Eq. (|8|)) versus the lattice coupling strength g for a 
two-dimensional CML with local map Eq. (^ on a square lattice of size L = 128 after an integration 
time of r = 3 X 10^ The magnetization bifurcates from a nonzero value at the critical value 
gc = 0.205. 

FIG. 4. Log-linear plot of the exponentially decaying two-point correlation function Eq. ^ 
of the two-dimensional CML Eq. (^ with local map Eq. (|^ on a square lattice for lattice sizes 
of L = 64, • ■ ■ , 1024. The coupling constant was g = 0.199 and the integration time was T = 3 x 10^. 
A reasonably accurate estimate of the correlation length ^2 — 6 is found only for sizes L > 256. The 
arrows indicate the range over which a least-squares linear fit was used to extract the length ^2- 
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FIG. 5. Plot of the two-point correlation function Eq. (^) for the magnitudes \yj\ of the field 
for the 2d periodic CML with local map Eq. (^) on a square lattice with g = 0.200. A system size 
L = 64 and an integration time of T = 5 x 10^ were used. The correlation length i^™'^^ ~ 1.4 is 
smaller than the length ~ 6 based on the values y*. 

FIG. 6. Log-linear plot showing the exponential decay of the mutual information func- 
tion I{Ax) as a function of the spatial separation Ax for the 2d periodic CML on a square lattice 
with local map Eq. coupling constant g = 0.204, system size L = 256, and effective integration 
time T = 128,000 (following a transient time of 100,000). The arrows indicate the range over 
which a linear least-squares fit was used to calculate the correlation length ^/ ~ 11.7. A simple 
binning procedure with 32 bins was adequate so that the more sophisticated algorithms of Eraser 
et al were not needed. 

EIG. 7. Two-point correlation length ,^2 versus the mutual information correlation length (^/ 
for a 2d periodic CML on a square lattice with local map Eq. (^). Lattices sizes of L = 1024 
and L = 256 and integration times of T = 300, 000 and T = 128, 000 were used respectively when 
calculating ^2 and The points correspond to the parameter range 0.1900 < g < 0.2045 with the 
smaller g values occurring to the left. The values of .^2 and were obtained from plots similar to 
Fig. I and Fig. | 

FIG. 8. (a) Divergence of ^2 for the 2d periodic CML on a square lattice with local map 
Eq. @. The lattice size was L = 1024 and the integration time was T = 3 x 10^ based on an en- 
semble average over 64 realizations. The correlation length diverges with critical exponent ~ —1.0 
at the extrapolated point g = 0.2054 which coincides within numerical accuracy with the pitchfork 
bifurcation of the magnetization, (b) The dimension correlation length ^5 (Eq. (|2|)) versus param- 
eter g for the same parameter values, based on dimension densities calculated from extensive chaos 
curves like Fig. ^. This length scale is roughly one lattice size, varies smoothly with g, and attains 
a maximum value for g ~ 0.200. 
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FIG. 9. Lyapunov exponents Aj versus the intensive index x = ijl? for the 2d periodic 
CML on a square lattice with local map Eq. ^ for the lattice size L = 24, for an integration 
time T = 50,000, and for lattice coupling values g = 0.19 and g = 0.22 which straddle the 
transition at g = 0.205. The exponents form an approximately continuous function that changes 
smoothly with parameter g. Only the first 512 of the 576 possible exponents are shown. 

FIG. 10. (a) The Lyapunov fractal dimension density (5, (b) the entropy density /i, and (c) the 
largest Lyapunov exponent Ai versus the coupling constant g in the vicinity of the ferromagnetic 
transition of gc = 0.205 according to the magnetization curve in Fig. |^ for the 2d CML on a square 
lattice with local map Eq. (^). System sizes of at least L = 16 were used with a fixed integration 
time T = 5 X 10^ and an ensemble average over p = 64 systems. The quantities 8 and h have local 
minima (although at slightly different g values) to the left of where the magnetization bifurcates. 
The error bars of each data point are about the size of the points themselves except for those of Ai 
which are much larger. 

FIG. 11. Magnetization M (Eq. dH)) versus coupling strength g for the 2d periodic CML on a 
square lattice but with the less chaotic lattice map given by Eq. (||). Other parameter values were 
the sizes, L = 64 and L = 128 and the total integration time, T = 3 x 10^. The bifurcation occurs 
at approximately gc ~ 0.168. 

FIG. 12. Lyapunov dimension density 6 and entropy density h versus coupling constant g 
for the 2d periodic CML with local map Eq. (|5|), for lattice size L = 32 and for an integra- 
tion time T = 10^. The largest Lyapunov exponent Ai = 0.095 remains constant in the range 
0.160 < g < 0.180 and so is not plotted. For this substantially less chaotic local map, extrema in 
the dynamical quantities do not occur. 
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FIG. 13. (a) Magnetization M (Eq. (P)) versus coupling strength g for 2d periodic CML on 
a hexagonal lattice with local map Eq. with integration time T = 1 x 10^ time units, for a 
system size L = 64, and for an ensemble average of p = 64 elements. A ferromagnetic transition 
occurs at g = gc ^ 0.12. (b) Dimension correlation length versus the lattice coupling constant g 
for the same parameters except with a system size L = 16 and integration time T = 5 x 10^. The 
maximum in S^g does not coincide with the bifurcation of the magnetization in (a). Over this range, 
the Lyapunov dimension decreases from a value of 106 to a minimum of D ~ 100 at g = 0.14. 

FIG. 14. (a) Magnetization M (Eq. (|8|)) versus lattice coupling strength g for the 3d periodic 
CML with local map Eq. (Q) on a cubic lattice. A system size of L = 16, a total integration 
time of T = 1 X 10^ time units, and an ensemble average of p = 64 realizations were used. The 
magnetization bifurcates to a nonzero value at g = gc ~ 0.11. (b) The dimension correlation 
length versus g for the same parameters, calculated on a lattice of size L = 6 and for a time of 
T = 1 X 10^ time units. The corresponding dimension D varies from 120 down to a minimum of 
about 105 at g = 0.15. 

FIG. 15. Variation of the correlation length the correlation length of magnitude fluctuations 
and the dimension correlation length with the radius of spatial coupling r for the Id periodic 
CML Eq. (^) with local map Eq. (Q). The system sizes ranged from L = 64 to L = 1024 with 
larger systems being used for larger radii. The integration time was T = 300, 000 for all systems. 
The length ^5 ~ 5.3 for r = 5 is substantially larger than the maximum value of ^5 ~ 1.2 found for 
the 2d Miller-Huse CML with nearest neighbor coupling. 
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Figure 1. O'Hern et al. Lyapunov Spectral Analysis.. 




Figure 2. O'Hern et al. Lyapunov Spectral Analysis 
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Figure 3. O'Hern et al. Lyapunov Spectral Analysis.. 




Figure 4. O'Hern et al. Lyapunov Spectral Analysis.. 
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Figure 5. O'Hern et al. Lyapunov Spectral Analysis.. 




Figure 6. O'Hern et al. Lyapunov Spectral Analysis.. 




Figure 7. O'Hern et al. Lyapunov Spectral Analysis.. 




Figure 8. O'Hern et al. Lyapunov Spectral Analysis.. 




Figure 9. O'Hern et al. Lyapunov Spectral Analysis.. 
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Figure 10. O'Hernetal. Lyapunov Spectral Analysis. 
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Figure 1 1 . O'Hern et al. Lyapunov Spectral Analysis.. 
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Figure 12. O'Hernetal. Lyapunov Spectral Analysis. 




Figure 13. O'Hernetal. Lyapunov Spectral Analysis.. 
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Figure 14. O'Hern et al. Lyapunov Spectral Analysis. 




Figure 15. O'Hernetal. Lyapunov Spectral Analysis.. 



